Loop dynamics and the evolution of enzyme activity.

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and ß-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Insights into the importance of WPD-loop sequence for activity and structure in protein tyrosine phosphatases.

Protein tyrosine phosphatases (PTPs) possess a conserved mobile catalytic loop, the WPD-loop, which brings an aspartic acid into the active site where it acts as an acid/base catalyst. Prior experimental and computational studies, focused on the human enzyme PTP1B and the PTP from Yersinia pestis, YopH, suggested that loop conformational dynamics are important in regulating both catalysis and evolvability. We have generated a chimeric protein in which the WPD-loop of YopH is transposed into PTP1B, and eight chimeras that systematically restored the loop sequence back to native PTP1B. Of these, four chimeras were soluble and were subjected to detailed biochemical and structural characterization, and a computational analysis of their WPD-loop dynamics. The chimeras maintain backbone structural integrity, with somewhat slower rates than either wild-type parent, and show differences in the pH dependency of catalysis, and changes in the effect of Mg2+. The chimeric proteins' WPD-loops differ significantly in their relative stability and rigidity. The time required for interconversion, coupled with electrostatic effects revealed by simulations, likely accounts for the activity differences between chimeras, and relative to the native enzymes. Our results further the understanding of connections between enzyme activity and the dynamics of catalytically important groups, particularly the effects of non-catalytic residues on key conformational equilibria.

Allosteric rescue of catalytically impaired ATP phosphoribosyltransferase variants links protein dynamics to active-site electrostatic preorganisation.

ATP phosphoribosyltransferase catalyses the first step of histidine biosynthesis and is controlled via a complex allosteric mechanism where the regulatory protein HisZ enhances catalysis by the catalytic protein HisGS while mediating allosteric inhibition by histidine. Activation by HisZ was proposed to position HisGS Arg56 to stabilise departure of the pyrophosphate leaving group. Here we report active-site mutants of HisGS with impaired reaction chemistry which can be allosterically restored by HisZ despite the HisZ:HisGS interface lying ~20 Å away from the active site. MD simulations indicate HisZ binding constrains the dynamics of HisGS to favour a preorganised active site where both Arg56 and Arg32 are poised to stabilise leaving-group departure in WT-HisGS. In the Arg56Ala-HisGS mutant, HisZ modulates Arg32 dynamics so that it can partially compensate for the absence of Arg56. These results illustrate how remote protein-protein interactions translate into catalytic resilience by restoring damaged electrostatic preorganisation at the active site.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (ßα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis.

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (µs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (ßα)8-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (ßα)8-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Exploiting enzyme evolution for computational protein design.

Recent years have seen an explosion of interest in understanding the physicochemical parameters that shape enzyme evolution, as well as substantial advances in computational enzyme design. This review discusses three areas where evolutionary information can be used as part of the design process: (i) using ancestral sequence reconstruction (ASR) to generate new starting points for enzyme design efforts; (ii) learning from how nature uses conformational dynamics in enzyme evolution to mimic this process in silico; and (iii) modular design of enzymes from smaller fragments, again mimicking the process by which nature appears to create new protein folds. Using showcase examples, we highlight the importance of incorporating evolutionary information to continue to push forward the boundaries of enzyme design studies.

The N-terminal Helix-Turn-Helix Motif of Transcription Factors MarA and Rob Drives DNA Recognition.

DNA-binding proteins play an important role in gene regulation and cellular function. The transcription factors MarA and Rob are two homologous members of the AraC/XylS family that regulate multidrug resistance. They share a common DNA-binding domain, and Rob possesses an additional C-terminal domain that permits binding of low-molecular weight effectors. Both proteins possess two helix-turn-helix (HTH) motifs capable of binding DNA; however, while MarA interacts with its promoter through both HTH-motifs, prior studies indicate that Rob binding to DNA via a single HTH-motif is sufficient for tight binding. In the present work, we perform microsecond time scale all-atom simulations of the binding of both transcription factors to different DNA sequences to understand the determinants of DNA recognition and binding. Our simulations characterize sequence-dependent changes in dynamical behavior upon DNA binding, showcasing the role of Arg40 of the N-terminal HTH-motif in allowing for specific tight binding. Finally, our simulations demonstrate that an acidic C-terminal loop of Rob can control the DNA binding mode, facilitating interconversion between the distinct DNA binding modes observed in MarA and Rob. In doing so, we provide detailed molecular insight into DNA binding and recognition by these proteins, which in turn is an important step toward the efficient design of antivirulence agents that target these proteins.

How abasic sites impact hole transfer dynamics in GC-rich DNA sequences.

Changes in DNA charge transfer properties upon the creation of apurinic and apyrimidinic sites have been used to monitor DNA repair processes, given that such lesions generally reduce charge transfer yields. However, because these lesions translate into distinct intra and extrahelical conformations depending on the nature of the unpaired base and its DNA context, it is unclear the actual impact of such diverse conformations on charge transfer. Here we combine classical molecular dynamics, quantum/molecular mechanics (QM/MM) calculations, and kinetic Monte Carlo simulations to investigate the impact of abasic sites on the structure and hole transfer (HT) properties of DNA. We consider both apurinic and apyrimidinic sites in polyG and polyGC sequences and find that most situations lead to intrahelical conformations where HT rates are significantly slowed down due to the energetic disorder induced by the abasic void. In contrast, the presence of an unpaired C flanked by C bases leads to an extrahelical conformation where stacking among G sites is reduced, leading to an attenuation of electronic couplings and a destabilization of hole states. Interestingly, this leads to an asymmetric HT behavior, given that the 5' to 3' transfer along the G strand is slowed down by one order of magnitude while the opposite 3' to 5' transfer remains similar to that estimated for the reference polyG sequence. Our simulations thus suggest that electrochemical monitoring of the DNA repair process following changes in charge transfer properties can miss repair events linked to abasic sites adopting extrahelical conformations.

Determination of the protonation preferences of bilin pigments in cryptophyte antenna complexes.

The light-harvesting mechanisms of cryptophyte antenna complexes have attracted considerable attention due to their ability to exhibit maximal photosynthetic activity under very low-light conditions and to display several colors, as well as the observation of vibronic coherent features in their two-dimensional electronic spectra. However, detailed investigations on the interplay between the protein environment and their light-harvesting properties are hampered by the uncertainty related to the protonation state of the underlying bilin pigments. Here we study the protonation preferences of four types of bilin pigments including 15,16-dihydrobiliverdin (DBV), phycoerythrobilin (PEB), phycocyanobilin (PCB) and mesobiliverdin (MBV), which are found in phycoerythrin PE545 and phycocyanin PC577, PC612, PC630 and PC645 complexes. We apply quantum chemical calculations coupled to continuum solvation calculations to predict the intrinsic acidity of bilins in aqueous solution, and then combine molecular dynamics simulations with empirical pKa estimates to investigate the impact of the local protein environment on the acidity of the pigments. We also report measurements of the absorption spectra of the five complexes in a wide range of pH in order to validate our simulations and investigate possible changes in the light harvesting properties of the complexes in the range of physiological pH found in the lumen (pH ∼ 5-7). The results suggest a pKa > 7 for DBV and MBV pigments in the α polypeptide chains of PE545 and PC630/PC645 complexes, which are not coordinated to a negatively charged amino acid. For the other PEB, DBV and PCB pigments, which interact with a Glu or Asp side chain, higher pKa values (pKa > 8) are estimated. Overall, the results support a preferential population of the fully protonated state for bilins in cryptophyte complexes under physiological conditions regardless of the specific type of pigment and local protein environment.

Single Amino Acid Mutation Controls Hole Transfer Dynamics in DNA-Methyltransferase HhaI Complexes.

Different mutagenic effects are generated by DNA oxidation that implies the formation of radical cation states (so-called holes) on purine nucleobases. The interaction of DNA with proteins may protect DNA from oxidative damage owing to hole transfer (HT) from the stack to aromatic amino acids. However, how protein binding affects HT dynamics in DNA is still poorly understood. Here, we report a computational study of HT in DNA complexes with methyltransferase HhaI with the aim of elucidating the molecular factors that explain why long-range DNA HT is inhibited when the glutamine residue inserted in the double helix is mutated into a tryptophan. We combine molecular dynamics, quantum chemistry, and kinetic Monte Carlo simulations and find that protein binding stabilizes the energies of the guanine radical cation states and significantly impacts the corresponding electronic couplings, thus determining the observed behavior, whereas the formation of a tryptophan radical leads to less efficient HT.

Molecular basis of the exciton-phonon interactions in the PE545 light-harvesting complex.

Long-lived quantum coherences observed in several photosynthetic pigment-protein complexes at low and at room temperatures have generated a heated debate over the impact that the coupling of electronic excitations to molecular vibrations of the relevant actors (pigments, proteins and solvents) has on the excitation energy transfer process. In this work, we use a combined MD and QM/MMPol strategy to investigate the exciton-phonon interactions in the PE545 light-harvesting complex by computing the spectral densities for each pigment and analyzing their consequences in the exciton dynamics. Insights into the origin of relevant peaks, as well as their differences among individual pigments, are provided by correlating them with normal modes obtained from a quasi-harmonic analysis of the motions sampled by the pigments in the complex. Our results indicate that both the protein and the solvent significantly modulate the intramolecular vibrations of the pigments thus playing an important role in promoting or suppressing certain exciton-phonon interactions. We also find that these low-frequency features are largely smoothed out when the spectral density is averaged over the complex, something difficult to avoid in experiments that underscores the need to combine theory and experiment to understand the origin of quantum coherence in photosynthetic light-harvesting.